Quantum state tomography of dissociating molecules 
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Using tomographic reconstruction we determine the complete internuclear quantum state, repre- 
sented by the Wigner function, of a dissociating I2 molecule based on femtosecond time resolved 
position and momentum distributions of the atomic fragments. The experimental data are recorded 
by timed ionization of the photofragments with an intense 20 fs laser pulse. Our reconstruction 
method, which relies on Jaynes' maximum entropy principle, will also be applicable to time resolved 
position or momentum data obtained with other experimental techniques. 



Quantum state tomography derives its name from the 
tomographic technique in medical diagnostics by which 
three-dimensional images of the inner parts of an object 
or a person can be derived from two-dimensional NMR or 
X-ray pictures obtained from different directions, a tech- 
nique for which Cormack and Hounsfield were awarded 
the Nobel prize in physiology and medicine in 1979. In 
quantum physics, we often deal with a phase space, char- 
acterizing the position and velocity distribution of par- 
ticles, and the aim of quantum state tomography is to 
determine this distribution from only position and mo- 
mentum measurements on the particles. As an example, 
the oscillatory motion of a particle in a quadratic po- 
tential is described as a simple rotation in phase space, 
and hence the mathematical Radon-transform technique 
0, used in physiology and medicine can be used to ex- 
tract the phase space distribution of such a particle if 
only position measurements are made sufficiently often 
over a time corresponding to the oscillation period. This 
technique has been succesfully applied to analyze bound 
molecular states, where the internuclear distance oscil- 
lates harmonically Q, and it has been used for non- 
classical states of light 0, 0, which are formally de- 
scribed as harmonically trapped particles. A number of 
extensions have been made to incorporate anharmonici- 
ties in the trapping potential, and to use mathematical 
reconstruction techniques which match specific detection 
schemes [U |H Hi ■ In the present work we determine the 
internuclear quantum state of dissociating J 2 molecules. 
The atomic fragments are created in a state with posi- 
tive energy, and they fly apart with no confining force 
between them. The situation thus differs from the one 
of trapped motion, and it calls for a new theoretical ap- 
proach. 

Experimentally, I2 molecules in a molecular beam are 
photodissociated by irradiation with a 100-fs-long pump 
pulse centered at 488 nm (see Fig. CJ). To dete rmmc 
the internuclear separation (position) and velocitiy (mo- 
mentum) of the dissociating molecules at time r after 
excitation the atomic / fragments are ionized by a 20- 
fs-long probe pulse, and the velocities of the I + ions are 
measured by a 2-dimensional ion detector 0] . If only one 



of the iodine atoms of a dissociating molecule is ionized 
the I + velocity is equal to half of the internuclear veloc- 
ity, approximately 4.0 A/ps. On the 2-D image, shown 
111 Fig. [U these ions constitute the innermost pair of half 
rings If, instead, both iodine atoms are ionized 

the velocity of the / + ions is increased due to their in- 
ternal Coulomb repulsion. On the 2-D image these ions 
constitute the outermost pair of half rings Us- 
ing Coulomb's law the / + velocity distribution from this 
channel gives directly the internuclear distribution [9j . 
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FIG. 1: Diagram of the experimental setup showing the 
molecular beam crossed at 90° by the pump and the probe 
laser beams. Both laser beams are polarized vertically. The 
I + ions produced are pushed by a weak electrostatic field (not 
shown) towards a 2-D ion detector, where their positions are 
used to determine the internuclear velocity distribution 
channel) and the internuclear separation distribution 
channel). The image shown on the detector is recorded for a 
delay of 5 ps between the pump and the probe pulse. 

The atomic density is sampled on a discrete position 
grid with N grid points, where N = 51 in the example 
shown in Fig. 2. For the data presented here, the recon- 
struction is based on position distributions recorded at r 
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= 2, 3, 4, 5 picoseconds (ps) and the momentum distribu- 
tion (which is identical at the four different delays). In 
quantum tomography on harmonic oscillator states the 
fast and slow velocity components are slushing back and 
forth, and the time dependent position distribution maps 
the phase space distribution from all angles. In free ex- 
pansion, slow particles will not be able to catch up on 
fast ones, and we have access to only sheared pictures 
of the distribution where the high velocity components 
are shifted towards larger displacements than the low ve- 
locity components. This implies that a range of angular 
viewing angles is missing, and we cannot apply the Radon 
transform to reconstruct the phase space distribution. 

Furthermore, the data-sets are noisy due to the finite 
counts in different position intervals, and our task is to 
establish the most reliable estimate of p in best possi- 
ble (rather than perfect) agreement with the observed 
data. Following Drobny and Buzek we apply Jaynes' 
maximum entropy principle which at the same time pro- 
vides the smallest deviation from the measured data 
and the largest possible von Neumann entropy 0, ^| : 
S(p) = —Tr(phip). Drobny and Buzek have used the 
maximum entropy principle for quantum state tomogra- 
phy on harmonically trapped atoms detected at too few 
instants of time to allow a normal reconstruction [13, 
and we have reformulated their method to deal with free 
atoms. 

To calculate the von Neumann entropy, one may use 
the familiar expression in terms of the eigenvalues pj of 
the density matrix, S = Pj^ n Pj-i but as it was 

shown by Jaynes [llj, it is not necessary to compute 
the entropy of p in order to use the maximum entropy 
principle: The density matrix with largest entropy which 
conforms best with measured values a, for a set of ob- 
servables Ai, can be written on the form 

P= ^exp(-^AiA), (1) 

i 

where Ai are variables that have to be adjusted to satisfy 
the agreement with the measured data, and Z normalizes 
the density matrix to unit trace. 

A position distribution is the set of expectation val- 
ues of projection operators on position eigenstates \x) (x\. 
We determine the position distribution at different times, 
and we need a formal representation of the projection 
operator \x,t)(x, t\ where \x,t) is a position eigenstate 
at time t. For free particles the momentum p = hk is 
a conserved quantity, and it is convenient to represent 
the projection operators in a momentum representation, 
since momentum eigenstates accumulate only a trivial 
phase factor with time, \k,t) = exp(ihk 2 t/2M)\k,t = 0): 

\x}(x\ = J dkidk 2 \ki) (ki\x) (x\k 2 ) (k 2 \ 

= — / dk 1 dk 2 exp(i(k 2 - k 1 )x)\k i )(k 2 \. (2) 
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Both the spatial distribution at various times and the mo- 
mentum distribution, which is independent of time, are 
measured, and the observables Ai in Eq.(l) are precisely 
the operators \x,t)(x,t\ and \k)(k\, which are now given 
cxplicity as N by N matrices in a discrete basis. The 
exponential in Eq.(l) is understood as a matrix power 
series, and it has to be evaluated numerically because the 
momentum projection operators and the position projec- 
tion operators at different times do not commute. 

In our numerical reconstruction we used four position 
distributions and one momentum distributions, and we 
thus identified 5-51 variables Ai so that the expectation 
values (Ai) =Tr(pAi) are as close as possible to the mea- 
sured distributions. We quantify the agreement with the 
measured data a; by the sum of the squares of all devia- 
tions A = J2i( a i ~ (Ai)) 2 - A numerical routine identifies 
the global minimum of A as a function of the Ai vari- 
ables, and we thus obtain the explicit expression for the 
density matrix in Eq.(l). 

To illustrate our reconstruction of the internuclear 
quantum state of the dissociated molecules we use the 
Wigner function W(x,p) ^3] rather than the N by N 
complex density matrix, discussed above. The Wigner 
function is obtained by a Fourier transform with respect 
to the difference between the two momentum arguments 
in the density matrix. While still fully characterizing 
the quantum state the Wigner function is a real func- 
tion, and it has the advantage of presenting at one glance 
the joint position and momentum distribution. Figure 2 
shows the Wigner function, W(x,p), at the earliest (2 
ps) and the latest time (5 ps) for which position distri- 
butions were recorded. Note that we present the momen- 
tum dependence in units of the velocity v = p/M of the 
atomic fragments to facilitate comparison of the two pan- 
els. The two Wigner functions in Fig. 2 a) and b) have 
the same velocity distributions, but they are sheared dif- 
ferently in phase space due to the free motion of the two 
iodine photofragments. On the side panels the marginal 
position and velocity distribution, obtained by integrat- 
ing W(x,p) over velocity and position, respectively, are 
compared with the experimental results (open circles). 
The agreement is excellent, and similar good agreement is 
found between the measured data and the reconstructed 
state at the intermediate times at 3 and 4 ps. 

We stress that the Wigner function contains more in- 
formation than what can be obtained directly from a 
single measured position and velocity distribution. For 
instance, both the measured position and velocity distri- 
bution at 2 ps consists of two peaks (side panels of Fig. 
2 (A)), but without further information it is not possible 
to determine if the molecules at the larger internuclear 
separation (the peak around 13.6 A ) are also the ones 
moving with the largest internuclear velocity (the peak 
around 4.6 A/ps). The tomographic reconstruction of 
the Wigner function makes use of our knowledge of the 
position distributions measured at 3, 4 and 5 ps, and, 
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FIG. 2: Three dimensional surface plots showing the phase 
space distribution function W(x,p) (Wigner function) at 
two different times after the pump pulse dissociates the Ii 
molecule (tau=2ps in part (A) and tau = 5 ps in part (B)). 
The marginal position and momentum distributions of the 
dissociated molecules are shown on the side panels: The full 
curves are the reconstructed distributions (integrals over the 
Wigner function with respect to momentum and position, re- 
spectively), and the open circles are the measured distribu- 
tions. To provide a more intuitive view the momentum vari- 
able is presented in velocity units. 



dissociation by the pump pulse |9j. 

It has been emphasized in the literature that the 
Wigner function is not a real probability distribution be- 
cause, due to complimentarity and Heisenbergs's uncer- 
tainty relation, it is not meaningful to assign precise val- 
ues for both the position and velocity of a particle. This 
has as its most striking consequence that states exist for 
which W(x,p) attains negative values, but these nega- 
tive values always occur in phase space regions with area 
less than Planck's constant Ti — h/2n. Dissociation of 
molecules by two phase-locked femtosecond laser pulses 
will produce photofragments in a coherent superposition 
of two localized wave packets 0, and we are currently 
working on improvement of the experiments to enable re- 
construction of the Wigner function with a sufficient res- 
olution to display the negativities expected in this case. 
Such an experiment constitutes a fs time resolved ana- 
logue to double-slit atom interferometer studies Q. 

Quantum state tomography has been established in 
quantum optics as a remarkable diagnostics tool that has 
provided clear illustrations of the experimental ability to 
control and manipulate selected states of fundamental 
quantum systems such as a single field mode and single 
ions and atoms 0, H, 0> H> E3 • Our experimentally re- 
constructed density matrix or Wigner function provides 
the complete information about the non-trivial quantum 
state of the internuclear motion of the dissociating I2 
molecules. Although our experimental method is lim- 
ited to small molecules, we believe that quantum state 
tomography will be applicable also to study chemical re- 
actions in large molecules through time resolved position 
data obtained, e.g., by the emerging femtosecond elec- 
tron and X-ray diffraction techniques 0, 0, |2(], l2l| . 
This will open for comparisons with theoretical calcu- 
lations at a much more detailed level than probability 
distributions of single observables, which have formed 
the basis of comparison so far in femtosecond time re- 
solved chemical reaction dynamics. Finally, we note that 
modern femtosecond laser technology enables controlled 
shaping of both the electronic and the atomic structure of 
molecules throu gh irradiation with sequences of tailored 
laser pulses 0,130. To fully exploit the capabilities 
of such quantum manipulation it is crucial to completely 
characterize the molecular quantum states formed. 

We acknowledge the support from the Carlsberg Foun- 
dation and The Danish Natural Science Council. 



indeed, the two separated peaks in W(x,p) at 2 ps show 
that the iodine atoms roughly 13.6 A apart have a ve- 
locity distribution centered around 4.6 A/ps whereas the 
iodine atoms apparoximately 12.4 A apart have a veloc- 
ity distribution centered at 4.0 A/ps. The physical origin 
of the two internuclear velocity components is the inco- 
herent population of v = and v = 1 vibrational levels 
of the electronic ground state of the I2 molecules prior to 
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